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ABSTRACT 



Context. Interferometers provide accurate diameter measurements of stars by analyzing both the continuum and the lines formed 
in photospheres and chromospheres. Tests of the geometrical extent of the chromospheres are therefore possible by comparing the 
\ estimated radius in the continuum of the photosphere and the estimated radii in chromospheric lines. 

Aims. We aim to constrain the geometrical extent of the chromosphere of non-binary K giant stars and detect any spatial structures in 
I the chromosphere. 

Methods. We performed observations with the CHARA interferometer and the VEGA beam combiner at optical wavelengths. We 
observed seven non-binary K giant stars (J3 and r] Cet, 6 Crt, p Boo, {S Oph, 109 Her and ( Cep). We measured the ratio of the radii 
of the photosphere to the chromosphere using the interferometric measurements in the H„ and the Call infrared triplet line cores. 
For J3 Cet, spectro-interferometric observations are compared to an non-local thermal equilibrium (NLTE) semi-empirical model 
atmosphere including a chromosphere. The NLTE computations provide line intensities and contribution functions that indicate the 
relative locations where the line cores are formed and can constrain the size of the limb-darkened disk of the stars with chromospheres. 
We measured the angular diameter of seven K giant stars and deduced their fundamental parameters: effective temperatures, radii, 
luminosities, and masses. We determined the geometrical extent of the chromosphere for four giant stars (/? and rj Cet, S Crt and p 
Boo). 

Results. The chromosphere extents obtained range between 16% to 47% of the stellar radius. The NLTE computations confirm that 
the Can/849 nm line core is deeper in the chromosphere of ft Cet than either of the Can/854 nm and Call/866 nm line cores. We 
present a modified version of a semi-empirical model atmosphere derived by fitting the Ca II triplet line cores of this star. In four of 
our targets, we also detect the signature of a differential signal showing the presence of asymmetries in the chromospheres. 
Conclusions. It is the first time that geometrical extents and structure in the chromospheres of non-binary K giant stars are determined 
by interferometry. These observations provide strong constrains on stellar atmosphere models. 



Key words. Techniques: interferometric - Stars: chromospheres - Stars: fundamental parameters - Stars: atmospheres - Radiative 
transfer 
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1 1- Introduction pulsation mechanisms, e.g. mass loss driven by accoustic waves. 

Th ese two groups of giant stars were proposed for the first time 



_ , , • , ■ by Linskv & Haisch ( 1979) from an analysis of ultraviolet spec- 
Chromospheres are the outer region of stars characterized by , 7 , . , . . . , .. ,• 

. . F r • tra. The group of coronal stars (giant stars hotter than K2) dis- 

a positive temperature gradient and a departure from radiative , . . ... . . . . , , 

JT b , • , play emission lines that must ongmnate in chromospheres, tran- 

equilibnum. These properties are caused by a heating mecha- ... . , , . ,. ™ c 

..,,,.,.„ , ... , sition regions, and by implication coronae. The group of non- 

msm in the low density part of stellar atmospheres, which results , . , ■ . . , , • • i- 

. J F , , . . F ' coronal stars (giant stars cooler than K2) display emission lines 

in temperatures ranging between the minimum of the tempera- c , . . ° , on nnn ts r.- i_ 

r, , ™ ™ formed at temperatures cooler than 20,000 K, which can ong- 

ture profile and ~ 20,000 K. The sources of this heating have . . , . , , ti.t-itt-i.j--j- i- 

, . , , • , , r ■ r mate only in chromospheres. The Linsky-Haisch dividing line 

not been tightly constrained but surface convection and mag- . e , ■ , . ■ " 

. _ , , 6 J , . , . ., , „ 6 separates the groups of coronal and non-coronal stars in the 

netic fields seem to be the main mechanisms responsible for TT „ /ttt ,n ■ • , . ■• , 

. . , r , ... Herzsprung-Russell (HR) diagram and was extensively studied 

heating the upper atmospheres of giant stars hotter than K2. in the 80 ' s (Simon et al . 1982; Ha isch 1987; Brown & Carpenter 

For cooler giants (V - R > 0.8 , these chromospheres look m4 Carpenter et al 1985) Even if this division of iant stars 

more extended relative to the stellar radius than tor hot giants . . , . . ,. ^ , , . . . ... 

, T , _ „. . , . . ... „ .j . into two groups seemed too simplistic, hybrid giant stars with 

(V — R < 0.8) and might be associated with stellar winds and 
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Table 1 : List of the K g iant stars observed with CHARA/VEGA. The magnitudes in the K band come from t he 2 MASS cata- 
log (|Skrutskie et ai1l2006l) . Atmospheric parameters (T e g, \ogg and [Fe/H]) come from iLuck&Challeneri (119951) and lMcWilliaml 
( 199Qb- Av is the interstellar absorption in the V band, p is the corrected Hipparcos parallax (Ivan Leeuwenll2007l) . 



HD 


Name 


Spectral Type 


P 


m v 


m K 


V -R 


Teff 


logg 


[Fe/H] 


A v 








(mas) 








(K) 








4128 


ySCet 


KOIII 


33.86 + 0.16 


2.02 


-0.20 


0.72 


4750 ± 100 


2.45 


0.13 





6805 


??Cet 


K1III 


26.32 ±0.14 


3.45 


0.86 


0.83 


4425 ± 100 


2.05 


0.04 





98430 


sen 


KOIII 


17.56 + 0.19 


3.56 


0.95 


0.83 


4500 ± 120 


2.59 


-0.48 


0.01 


127665 


p Boo 


K3m 


20.37 +0.18 


3.59 


0.59 


0.92 


4260 ± 120 


2.22 


-0.17 





161096 


/3 0ph 


K2III 


39.85 + 0.17 


2.77 


0.27 


0.82 


4475 ± 100 


1.70 


0.00 


0.05 


169414 


109 Her 


K2III 


27.42 ± 0.40 


3.84 


1.09 


0.85 


4450 ± 120 


2.67 


-0.16 


0.04 


216228 


( Cep 


KOIII 


28.29 + 0.10 


3.53 


1.11 


0.83 


4770 ± 120 


2.97 


-0.12 






mixed coronae and signific ant wind activity have been revealed 
(lReimerslll98l lHalll (20081 and lAvresll2010l) . For example, the 
K3III giant star, 6 And, shows an unexpected presence of C 
IV in emission, which implies that it contains hot material 
about 100,000 K) , and evidence of a strong, high-velocity wind 
Judge et al] |1987l) . This hybrid state appears to represent the 
transition between coronal and non-coronal groups. 
Whatever the classification into either coronal, non-coronal, or 
hybrid groups, no real agreement has been reached about the ge- 
ometrical ext ent of the chromosphere of cool giant stars. Finally, 
Judge] (1 19871) argued the chromosphere extents do not change 
dramatically as a star crosses the Linsky-Haisch dividing line. 
Th is resul t is also s upported by model atmospheres developed 
by ICuntzl (Il990albl) . which predict their extent to be betwee n 
5% and 50% of the stellar radius (see Table 3 of lCuntzl lT990a). 
These model atmospheres are constructed using spectroscopic 
observations of C n, Mg n, and Ca n lines in the UV and IR wave- 
length ranges, forcing the radiative transfer of these lines to fit 
their cores formed in the chromosphere, transition, or coronae 
regions. 

It is apparent that the extents of the chromosphere and its 
physical conditions are still not fully understood for cool giant 
stars. The geometrical constraints on observed chromospheres 
are non-existant except for eclipsing binaries for which the ra- 
dius can be m easured in a spectral line as for example for the 
£ Aur system dEaton|[l993l) . It is therefore worthwhile to do in- 
terferometric observations of red giant stars in the cores of lines 
formed in chromospheres and to compare the resulting measured 
radii with those obtained from the continuum during the same 
run of observations. For this purpose we need an interferome- 
ter working in the v isible or near-IR such as CHARA/VEGA 
(Mou rardetal.1 l2QQ9h allowing us to observe the H a line, the 
Ca n triplet lines (849, 855 and 860 nm), and the nearby contin- 
uum. 



This work compares the measured radii in different spectral 
lines and the continuum for seven K giant stars, one of these 
stars being classified as a coronal star (J3 Cet). The observa- 
tions and data processing are described in Sect. 2. In Sect. 3, 
we present our estimates of the fundamental parameters of the 
program stars, followed in Sect. 4 by our study of the chromo- 
sphere of four giant stars. In Sect. 5, a semi-empirical model of 
the chromosphere of /3 Cet is used to interpret spectroscopic and 
interferometric data. Conclusions are presented in the last sec- 
tion. 



2. Interferometric observations 

2.1. Data 

Seven K giant stars, listed in Table Q] were observed at medium 
and high spectral res olution with the Visible spEctroGraph And 
polarimeter (VEGA, iMourard et al.ll2009h integrated within the 
CHARA array at Mount Wi lson Observatory (California, USA, 
ten Brummelaar et al. 1 12001 . We used two criteria to select the 
giant stars for the program: 

- the visibility in the stellar continuum at 790 nm had to be 
greater than 30%; 

- the stars had to be bright enough, i.e. my < 4 (in high spectral 
resolution mode). 

These two criteria ensured the optimal operation of VEGA when 
no fringe tracker is used. 

Apart from [5 Ophicus, i Cephei, and 109 Herculis, which 
were observed in the continuum only, all stars were observed 
in both the continuum and chromospheric spectral lines. We 
selected four lines with chromospheric cores: the H a Balmer 
line and the Can infrared tripet lines (849.8 nm, 854.2 nm and 
866.2 nm). Details of the observations can be found in Table|2] 
The shortest baseline of the CHARA array (S1S2, 33 m) was 
used during the observations and the data were recorded in two 
spectral bands of the continuum simultaneously using the two 
detectors of VEGA (blue and red detectors). Calibrator stars 
were also observed in the continuum in order to calibrate the 
measurements. We used the following sequence of observations 
cal-target-cal-target-ccil when only one calibrator was avail- 
able and the sequence call -cal2-target-cal2-call -target-call - 
cal2 when two calibrators were available. We used the medium 
spectral resolution mode of VEGA (R = 5, 000) for the observa- 
tions in the continuum. 

The high spectral resolution mode of VEGA (R = 30, 000) 
was used for the observations of spectral lines for which calibra- 
tors were not required because the visibility in the spectral lines 
were calibrated by the measurements in the continuum, close to 
the spectral lines (see Sect. 12.21 1. 

We selecte t he calibrators using t he SearchCal tooQ devel- 
oped at JMMC dBonneau et al.ll2006h . providing an estimate of 
the limb-darkened (LD) angular diameter (0u>). The uniform- 
disk (UD) angular diameter (Sj/d) is required to estimate the 
transfer function of the instrument at each wavelength. The 
UD angular diameter of each calibrator i s derived using 6j,n 
with the linear LD coefficients given by IClaret et all {1995 ) 
and iDiaz-Cordoves et aD (1 19951) . For each calibrator, Table [3] 
presents 9ud at 620 nm and 790 nm. 



1 available at: http://www.jmmc.fr/searchcal_page.htm 
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Table 2: Journal of the observations. Aq is the central wavelength 
of the recorded band. The last column presents the spectral band- 
width A/l used in the data processing (see Sect. 12.2b . 



Name 


Type 


Epoch 


^0 

(nm) 


HA 


/J Ceti 


target 


10/08/26 


790.0 


10 nm 






10/08/26 


620.0 


15 nm 






10/07/31 


Can 849.8 


0.25 A 






10/07/^1 


Ca TT R54 ? 


0.25 A 






10/07 HI 

1U/U/ / Jl 


Pa n RfSfS ? 


0.25 A 




cat. 


1 n /ns /if,. 

lU/Uo/ ZO 


/yu.u 


1 n nm 

iu nm 






1U/ UO/ 






r] Ceti 


target 


10/09/16 


790.0 


10 nm 






10/09/16 


620.0 


15 nm 






10/00/1 7 

J-U/ U7/ 1 / 


Ca TT R4Q R 
V_-d. 11 0^".0 


1 A 






J-U/ U7/ 1 U 


Pa tt ? 


1 A 






1 A /no / 1 "7 

1U/U9/ 1 / 


L.ail ooo.z 


1 A 

1 A 




cal. 


1 n /no / 1 a 


/yu.u 


10 nm 






1 n /no / 1 & 


ozu.u 


15 nm 


TT* ■ 

o Cratens 


target 


IU/Uj/Uj 


/yu.u 


10 nm 






lU/UO/Uj 


ozu.u 


15 nm 






10/05/05 


/~i or /i ^ 

Can 854.2 


l A 






10/05/05 


H„ 656.2 


lA 


HD100889 


cal. 


10/05/05 


790.0 


10 nm 






10/05/05 


620.0 


15 nm 


p Bootis 


target 


10/06/26 


790.0 


10 nm 






10/06/26 


620.0 


15 nm 






10/07/31 


Can 854.2 


1 A 


HD 143 894 


cal. 


10/06/26 


790.0 


10 nm 






10/06/26 


620.0 


15 nm 


HD108382 


cal. 


10/06/26 


790.0 


10 nm 






10/06/26 


620.0 


15 nm 


P Ophicus 


target 


10/09/16 


790.0 


10 nm 






10/09/16 


620.0 


15 nm 


HD152614 


cal. 


10/09/16 


790.0 


10 nm 






10/09/16 


620.0 


15 nm 


109 Herculis 


target 


10/07/31 


790.0 


10 nm 






10/07/31 


620.0 


15 nm 


HD168151 


cal. 


10/07/31 


790.0 


10 nm 






10/07/31 


620.0 


15 nm 


HD 166230 


cal. 


10/07/31 


790.0 


10 nm 






10/07/31 


620.0 


15 nm 


i Cephei 


target 


10/11/10 


790.0 


10 nm 






10/11/10 


620.0 


15 nm 


HD3360 


cal. 


10/11/10 


790.0 


10 nm 






10/11/10 


620.0 


15 nm 



2.2. Data processing 

The data processing of VEGA is composed of two parts. First, 
the data in the stellar continuum (at medium spectral resolution) 
are processed with the power spectral method giving the squared 
visibilities. Second, the processing of the data in spectral lines is 
based on the cross-spectrum method, which provides differential 
visibilities and phases across the lines (see Mo urard et al.1 2009 
for details). 

For this analysis and in the case of observations at medium 
spectral resolution, we divide the whole spectral band recorded 
by the red detector into four spectral channels of 10 nm centered 
on the wavelengths 775 nm, 785 nm, 795 nm, and 805 nm. For 
the blue detector, the continuum is visible only at 625.5 nm, we 
then use only one spectral channel (with a spectral bandwidth 
of 15 nm). The processing of these five spectral bands give the 
squared visibilities used to constrain the LD angular diameter 
(see Sect. Em 



Table 3: Limb-darkened On, and uniform-disk Oud angular di- 
ameters of calibrators. 



Name G LD 8 UD at 620 nm 9 UD at 790 nm 



(mas) (mas) (mas) 

HD3360 0.290 ± 0.020 0283 0284 

HD166230 0.420 + 0.029 0.393 0.398 

HD168151 0.679 + 0.047 0.663 0.670 

HD152614 0.312 ±0.022 0.301 0.303 

HD108382 0.394 ±0.027 0.376 0.380 

HD143894 0.379 ± 0.026 0.362 0.365 

HD100889 0.305 ±0.021 0.293 0.296 

HD222345 0.511 ±0.035 0.485 0.491 

HD225132 0.351 ±0.024 0.338 0.341 



The cross-spectrum method is applied to the data recorded 
in the high spectral resolution mode. We compute the complex 
differential visibility between a large spectral channel used as 
reference (centered at A\) and a narrow spectral band (centered 
at A2) sliding in the reference spectral channel 

Vdiff = V Al V M exp [-i(<p M - <p A2 )], (1) 

where Va i and 0^. represent the visibility and the phase of the 
fringe patterns at the wavelength Aj. The width of the sliding nar- 
row spectral band (A/^ hereafter) is chosen to ensure a sufficient 
signal-to-noise ratio (S /N) for the estimation of the complex dif- 
ferential visibility. Therefore, it depends on the stellar magni- 
tude. We use the following spectral bandwidths: AA2 = 0.25 A 
for p Cet and AA 2 - 1 A for 6 Crt, p Boo, and rj Cet. 
The cross-spectrum method recovers both the modulus of Vdiff 
(the differential visibility) and its argument (the differential 
phase A</>i2 = <px x - <Px^)- As no calibrator has been observed 
for the high spectral resolution mode, we developed a special- 
ized processing method to calibrate the visibility in the narrow 
spectral band. The sequence of processing removing the instru- 
mental/atmospheric signature was: 

1 . Estimation of Vdiff using the cross-spectrum method. 

2. Fitting of V^m wit h a simple model as described in 
Mourard et al. (2009). This model describes the instrumen- 
tal/atmospheric signature. 

3. Normalization of Vdiff with the fitted model. 

We then obtained V^, normalized to one in the continuum. 

4. Using the angular diameter deduced from our observations 
in the continuum, we computed the theoretical visibility to 
rescale Va 2 (see Fig. [T}. 



3. Fundamental parameters 

3.1. Angular diameter 

The squared visibilities in all spectral channels are fitted with 
a model of LP di sk, w hich is constructed u s ing th e laws of 
IClaret et alJ dl995l) and iDiaz-Cordoves et all dl995h . Table H 
presents the LD angular diameters Bid- Comparison with a 
previously published diameter is only possible for f3 Cet for 
which the diameter was estimated with the instrument VINCI 
at VLTI in the K band dRichichiet al.ll2009h . Our estimate for 
/? Cet (5.288 + 0.075 mas) agrees well with the VINCI one 
(5.329 + 0.005 mas). 
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Table 4: Fundamental parameters of the observed K giant stars. 



HD 


Name 


old 


/bol 




L/L e 


R/R Q 






(mas) 


erg ra s cm ) 


(K-) 








4128 


jSCeti 


5.288 ± 0.075 


5.10 


4838 ±70 


139.1 ± 


7.0 


16.78 ± 0.25 


6805 


rj Ceti 


3.698 ±0.160 


1.64 


4356 ± 55 


74.0 ± 


3.7 


15.10 ±0.10 


98430 


6 Crateris 


3.667 ± 0.022 


1.66 


4408 ± 57 


171.4 ± 


9.0 


22.44 ± 0.28 


127665 


p Bootis 


4.090 ±0.031 


1.76 


4298 ± 56 


131.9 ± 


6.8 


21.57 ±0.25 


161096 


P Ophicus 


4.606 ± 0.045 


3.04 


4621 ±62 


63.4 ± 


3.2 


12.42 ±0.13 


169414 


109 Herculis 


3.223 ± 0.034 


1.17 


4334 ± 59 


50.7 ± 


2.7 


12.63 ± 0.22 


216228 


l Cephei 


2.646 ± 0.048 


1.28 


4831 ±74 


49.6 ± 


2.5 


10.05 ±0.18 




STEP 1 and 2 



E °-« 

O 0.2 



OJ 0.8 
2 0.6 



STEP 4 



854.0 854.5 855.0 

Wavelength {ntn) 



Fig. 1: Illustration of the four steps of the processing of observa- 
tions at high spectral resolution in an absorption line. These data 
correspond to 77 Cet observed in the Ca n/854 nm line. Top: raw 
differential visibility with the simple model. Middle: normalized 
visibility in the sliding narrow spectral band. Bottom: calibrated 
visibility in the sliding narrow spectral band. 
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- H 


1 








1,1, 





4400 4600 
spectroscopic Tell (K) 



Fig. 2: Comparison of spectroscopic and interferometric estima- 
tions of T e ff for our program stars. The spectroscopic T e g come 
from Table Q] 



3.2. Effective temperature 

The effective temperature r e ff of a star, its LD angular diameter 
Old, and the bolometric flux /i,oi are related via the equation 



4/t 



bol 



eff 



Baines et al 
Di Benedetto et al 
Dyck et al 
Van Belle et al 
Mozurkewich et al 
this work 




Fig. 3: Comparison of our estimates of T e ff for our program stars 
with the temperature scale deduced from other interferometric 
measurements. 



where cr is the Stefan Boltzmann constant. The bolometric fluxes 
of each star are given in Table [4] and were computed us i ng the 



(2000) 



Chen et al. 



bolometric correction in K band from Houdashe lt et al. 
and th e interstellar absorption Ay from the work of Q 
(Il998h . For all the program stars, it turns out that the absorp- 
tion is very small because they are not in the direction of the 
Galactic plane and are closer than 100 pc. The interferometric 
r e ff obtained are presented in Table [4] We assume an error of 
5% in /b i, which is associated with the error in the photome- 
try. Therefore, the errors in the range from 55 K to 74 K, 
implying errors of about 1.5% in the interferometric T e $. 

A comparison of our interferometric T ^f estimates with 
the sp ectr oscopic r e g found in the literature dLuck & Challenerl 
1995 and McWilliaml l 19901) is presented in Fig. [2] showing that 
our interferometric measurements are within the spectroscopic 
error bars. In Fig. [3] we show that our T e ff estimates agree 
very well with the temperature scale r e ff versus (V - K) de- 
duced from the in terferometric r e ff e s timated by several authors 
for K giant sta rs (Baine s et al.ll2010|: | M ozurkewich et al. 2003; 
van B elle et al.|[l999t iDvck et all 1 19961: Idi Benedetto & R abbia 
1987). This figure shows that our program stars follow the fitted 
linear-regression relation and that ou r estimates are wi thin the 
spread of AT^ * 145 K proposed bv lDvck eTall (119961) . 

3.3. Position in the HR diagram 

The radius is obtained with the formula 



(2) R/R Q = 



0. 



LD 



9.305/? 



(3) 
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0.5 
2.5 



2.0 



1 5 




0.5 



Ceti 
^7] Ceti 

^.p Bootis 

^(3 Ophicus „ 
Cephei 
109 Herculia. 



- [Fe/H]=-0.35 




6 Crateris _ 



6 Crateris 



3.70 3.65 
loq(Teff) 



Fig. 4: HR diagram of the program stars and comparison with 
evolutionnary track models for different masses. Top: [Fe/H] = 
0.06; bottom: [Fe/H] = -0.35. 

Table 5: Estimated masses of the seven observed K giant stars 
deduced from the BaSTI evolutionnary tracks. 



Name 


Mass 




(Af E ) 



/3Ceti 


3.0- 


-3.5 


j] Ceti 


1.0- 


- 1.4 


6 Crateris 


1.0- 


- 1.4 


p Bootis 


1.0- 


- 1.4 


P Ophicus 


1.4- 


- 1.8 


109 Herculis 


0.8- 


- 1.0 


i Cephei 


1.8- 


-2.2 



where p is the paral lax in arcseconds co ming from the corrected 
Hipparcos catalog (Ivan Leeuwerfe 007) and Old is in mas. The 
luminosity is obtained from the absolute bolometric magnitude 
computed with the apparent magnitude in V band (see Table []]), 
the c orrected Hipparco s paral lax, and the bolometric correction 
from iHoudashelt et al] ([2000). Results for L/L Q and R/R Q are 
given in Table [4] Plotting these results in a HR diagram to- 
gether with evolutionary tracks results in a rough estimate of 
the mass of the program stars (see Fig. 0). We used evolution- 
ary track models at solar metallicity except for 6 Crt for which 
we use models at [Fe/H]=-0.35, as shown in the bottom panel of 
Fig. [4] The evolutionar y tracks come from the BaSTI database 
( Pietrinfer ni et alJ20 04). Table[5]presents the mass estimates de- 
duced from this comparison. For/? Cet, we deter mine a more ac- 
curate mass of 3.1M G using the CESAM code dMorelll 19971) . in 
good agreement with the results presented in Table [5] 



0.8 
>, 0.6 



^ 0.4 
> 



0.0 
30 
20 
10 



■a 



g-20 
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655.5 656.0 656.5 
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Fig. 5: Plots of the visibility and the differential phase versus 
wavelength in the case of observations around the H a line (dotted 
vertical line) for 5 Crateris. 



4. Chromosphere parameters 

4.1. Geometrical extent 

Visibility measurements in the core of the chromospheric lines 
give us direct estimates of the extent of the line-forming regions. 
For this purpose, we assume that the line cores are optically thick 
in the region where they are formed, i.e. the chromosphere. This 
is demonstrated in Sect. [5] for ft Cet for which we have a semi- 
empirical model of the extended atmo sphere. This hypothesis 
agrees with similar studies for the Sun dWhitdll962l) . Hence, in 
a first approximation, we can consider a UD model to derive the 
extent of the chromosphere. 

Applying the cross-spectrum method to reduce the data recorded 
in the high spectral resolution mode (see Sect. l2.2K we derive the 
visibility curves as presented in Figs. [5] [6] El and[8] 
For each giant star, we then deduce a UD angular diameter from 
the value of the visibility in the core of the lines. The results are 
presented in Table [6] which also gives the extent in linear stellar 
radius equal to the ratio of the LD angular diameter in the 
lines to the stellar LD angular diameter given in Table [4] In a 
first approximation for computing Old in the core of the lines, 
we use the same correcting factor OldI^ud ~ 1-05 given by the 
limb-darkening laws derived with the semi-empirical model of 
(5 Cet (see Sect©. 

In our program stars, only /? Cet belongs to the group of coronal 
stars, which means it is located on the left of th e Linsky-Haisch 
dividi ng line of the HR diagram (V — R = 0.72). ICarpenter et al.l 
(1 19851) also give the same classification for [} Cet. We find that 
the extents of the emitting regions of the Can triplet and H a 
lines range from 16% to 47% of the stellar radius. These results 
show that the extent of the chromospheric line-forming regions 
is similar for the four cool giant star s regardless of their classi- 
fication. This conclusi on agrees with Judge ( 1987), but not with 
ICarpenter et al.l (1 19851) who proposed the idea that coronal stars 
appear to have a very thin chromosphere of no more than 0. 1 % of 
the stellar radius. Moreover, our estimated chromosphere sizes 
agree well with the theoretical predictions of ICuntzl (11990a). 

4.2. Structures at the surface 

Information related to the position of the photocenter of the 
chromosphere can be deduced by measuring the phase of the 
visibility (also called differential phase) in the narrow spectral 
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Table 6: Extent of the chromospheric line-forming regions. The extents are given in UD angular diameter and in stellar radius. 



Name 


UD Angular diameter (6 UD in mas) 




Linear radius (in 






Ca n/849 nm 


Ca n/854 nm 


Ca n/866 nm H„ 


Can/849 nm 


Ca n/854 nm 


Ca n/866 nm 


H„ 


ySCeti 


5.82 ± 0.20 


6.40 ± 0.23 


6.58 ± 0.24 


1.16 + 0.04 


1.27 ± 0.05 


1.31+0.05 




7/ Ceti 


4.40 ±0.19 


4.87 ± 0.20 


5.05 ± 0.24 


1.25 + 0.05 


1.38 ± 0.06 


1.43 + 0.07 




6 Crateris 




5.10 + 0.30 


5.13 + 0.18 




1.46 ± 0.09 




1.47 + 0.05 


p Bootis 




4.86 + 0.28 






1.25 + 0.08 







bands. This observable has been intensively used to study the 
kine matic of circumstella r environments (e.g. iBerio et al.lll999l 
and iMeilland etaHl2007b . In Figs. [5] EJ and |8] we plot the 
differential phase with respect to the wavelength used for the 
observation of the chromospheric lines. Clear jumps in the dif- 
ferential phase at the center of the Ca n/854 nm line are revealed 
for three of the four program stars observed at high spectral res- 
olution. Only rj Cet does not exhibit an unambiguous signal. For 
Can/849 nm and H a lines, no jump is detected in the differential 
phase. For the Ca n/866 nm line, jumps seem to be present but 
their amplitude with respect to the phase precision is too low for 
a significant detection. 

The phase jumps can be explained either in terms of an asym- 
metric chromosphere or an asymmetric photosphere, which im- 
plies structure in the chromosphere and/or in the photosphere. 
Our observations provide phase measurements for only one 
baseline for each line preventing us from deriving the exact pho- 
tocenter positions. However, the fact that phase signatures are 
not present with the same amplitude in all lines may indicate that 
structures are present in the chromosphere and not only in the 
photosphere. The star 6 Crt perfectly illustrates this property be- 
cause an unambiguous phase jump (amplitude greater than 20° 
seen in the bottom left panel of Fig. [6]l is present in the core of 
the Ca n/854 nm line whereas the phase remains constant over 
the H a line (see Fig. [5]). 

Until now, only spectroscopic methods have been used to 
detect and study the presence of stru ctures in the chromo sphere 
of single K giant stars. For instance. iMontes et alj ((2000) stud- 
ied in detail the chromosphere, discriminating between different 
structures (plages, prominences, flares, and microflares), by an- 
alyzing the ratio of excess emission equivalent width (EW) of 
two Can infrared lines (£ , 8542/£'8498) or the rat i° of EWs of two 
B aimer lines (E^/E^). Our study shows that additional infor- 
mations could be provided by interferometry in the future to help 
us understand the structures in the chromosphere of single K gi- 
ant stars. In reality, the structure of the chromospheric network 
delineated by bundles of magnetic field lines could be spatially 
studied in the future. Combinations of interferometric observa- 
tions in several spectral lines will allow us to determine the char- 
acteristics (size, position, intensity) of the chromospheric struc- 
tures such as plages (part of the chromospheric network of bright 
emission associated with concentrations of magnetic fields) and 
prominences or filaments (dense clouds of material suspended 
above the stellar surface by loops of magnetic field). These per- 
spectives require interferometric observations with a more com- 
plete (m, v) coverage (several baselines with different lengths and 
orientations) than that available from these first observations. 



5. Semi-empirical model of /? Cet 

The star f3 Cet is the only one for which a model atmosphere 
including a chr omosphere and a transition region has been fitted 
(Erik sson et al.|[l983l) . The adopted atmospheric parameters for 
the target are given in Table [7] 



Table 7: Comparison of the model atmosphere parameters for (3 
Cet. Second column is the model geometry: plane-parallel (pp) 
or spherical (sph). Last column is the plateau temperature of the 
chromosphere T p . 



Model atmosphere 


geo. 


Teg 


l°gg 


[Fe/H] 


T P 






(K) 






(K) 


Eriksson et al. (1983) 


PP 


4900 


2.90 





5500 


This work 


sph 


4830 


2.45 





4500 



We note that lEriksson et al.l (Q983) used a pre-Hipparcos par- 
alax estimate of p — 61 mas , which is about twice the Hipparcos 
measurement (p = 34 mas. iPerrvman et afll 19971) . The star was 
thought to be closer, thus smaller and con sequently leading to a n 
overestimation of the surface gravity by lEriksson et al.l (|J983). 
We adopt a surface gravity based on our previously estimated 
interferometric radius (R = 16.8 R e ) and the mass deduced from 
the evolutionary tracks (M = 3 M Q ). The uncertainties in these 
two parameters give a surface gravity ranging from 2.45 to 2.55. 
We adopted a value of log e = 2.45, which is consis tent with the 
values determined by both Lu ck & Challenerl(ll995l) and lKovacsl 
(119831) from a spectroscopic analysis based on the ionization 
equilibri um of iron. We interp olate the MARCS0 model atmo- 
spheres dGustafsson et al.1120081) with spherical geometry for the 
atmospherical parameters T e $ = 4830 K, logg = 2.45, [Fe/H] 
= 0, and £ = 2 km s , all obtained with the interpolation code 
of T. Masseroifl For the interpolation, we note that we use four 
spherical model atmospheres corresponding to a stellar mass of 
2 M available in the MARCS grids. The theoretical MARCS 
model atmospheres do not include either chromosphere or tran- 
sition region. To re produce the core of th e Can triplet lines, 
we use the model of Eriksson et al. ( 1983) replacing its photo- 
spheric part by the interpolated MARCS model, hereafter called 
the hybrid model. 

The flux and intensity profiles are computed using the 
non-local t hermal eq uilibrium (NLTE) radiative transfer code 
M ULTI2.2 dCarlssonl 19861) with the Ca n model atom described 
in M erle et al.l(l201ll) and a model atmosphere corresponding to 
(3 Cet. We notice that under the assumption of local thermal equi- 
librium (LTE), the Ca n triplet lines appear with a core in emis- 
sion if we use a model atmosphere with a chromosphere. This 
justifies the use of NLTE radiative transfer in correctly treating 
the line formation. The MULTI2.2 code provides the contribu- 
tion functions of the flux showing where the different parts of 
the line are formed in the atmosphere. We see in Fig.[TT]that the 
continuum is formed in the deep photosphere and the cores of 
the Ca ii triplet lines between log m - -2 and -4 in the atmo- 
sphere, corresponding to the mean chromosphere (see top panel 
of Fig. |5J. We note that the contrib ution function in the contin- 
uum for the lEriksson et al.1 (l983) model is truncated. This is 
because their model does not consider sufficient depths of the 



2 available at: http://marcs.astro.uu.se 

3 available at: http://marcs.astro.uu.se/software.php 
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Fig. 6: Plots of the visibility and the differential phase versus wavelength in the case of observations around the Can/854 nm line 
(dotted vertical line) for four K giant stars. 
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Fig. 7: Plots of the visibility and the differential phase versus wavelength in the case of observations around the Can/849 nm line 
(dotted vertical line) for two K giant stars. 



photosphere, as shown on the top panel of Fig. [9] The maximum 
value of the contribution functions in the line cores of the Ca n 
triplet corresponds to the mean chromosphere characterized by 
a constant temperature (see Fig. [9j, i-e. the plateau temperature 
T p . The plateau tem perature was a parameter of the model of 
Eriks son etail d 19831) set to fit the Mg n k line flux from the IUE 
observations. From Fig.QTl we conclude that the chromosphere 
is optically thick in the core of the Ca n triplet lines. Using the in- 
terpolated MARCS model without chromosphere, the line core 
contribution functions are zero, confirming this conclusion. In 
addition, we note that in Fig.[TT] the Can/849 nm line is formed 
more deeply in the chromosphere than the other two lines. This 
prediction is clearly supported by our interferometric measure- 
ments (see Table[6]). We emphasize that we measure similar rel- 
ative chromosphere extents in the three components of the Ca n 



triplet for the only other star, rj Cet, for which we have interfer- 
ometric observations in the core of the Ca n triplet. 

We use MULTI2.2 to compute the LD laws presented in 
Fig. [TO] for the core of the Can triplet lines and for the close 
continuum. For this purpose, we modified the line transfer in 
MULTI2.2 code in order to compute 32 directions (fj. = cos 6) 
between and 1. Our result validates the approximation of 
the UD model for deriving the chromospheric radius, shown in 
Table[6] The computations invalidate the possibility that the line 
cores could be brighter in the limb than in the center-disk, as we 
thought at the beginning of the study. The shape of the LD law 
in the chromospheric line core is simil ar to the sh ape of the H a 
Balmer lines of the Sun derived by White ( 119621) . We then use 
our model to derive Old ■ For that, we fit to each spectral line a 
third-order polynomial law in fi to the intensity profile presented 
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Fig. 8: Plots of the visibility and the differential phase versus wavelength in the case of observations around the Ca n/866 nm (dotted 
vertical line) line for two K giant stars. 



in Fig. [10] These coefficients are then used to derive 610 by fit- 
ting a third-order LD law to our visibility measurements. We find 
that the radius of each Ca 11 line core is larger by the same fac- 
tor of ~ 5% (6> LD (849 nm) = 6.11 + 0.21 mas, 6» LD (854 nm) = 
6.72 + 0.24 mas, and 6> L d(866 nm) = 6.91 ± 0.25 mas) than the 
9ud presented in Table [6] 

Finally, we compute NLT E Ca 11 line profiles for three model 
atmospheres: (1) the model of Eri ksson et al.l (I 1 9831) with the in- 
correct surface gravity, (2) the interpolated MARCS model with- 
out chromosphere, and (3) the hybrid model. The NLTE flux 
profiles of Ca 11 triplet lines are compared with the observations 
of (3 Cet from the NARVAL instrument at TBlfl (R ~ 65, 000; 
S IN ~ 500). The project ed rotational velocity o f j3 Cet is es- 
timated to be ~ 5km s -1 (iMassarotti et al.l [2008). We included 
the convolution of the intrinsic flux profile with the rotation pro- 
file befor e degrading to t he sp ectrograph resolution. Using the 
model of Erik sson et al.l (Il983l) . we find evidence of emission 
near the core in Fig. [T_2] especially for the Ca 11/866 nm line. 
We interpret this emission as an indication of a higher temper- 
ature T p . As an exercice, we scale the chromos phere and tran- 
sition part (as proposed by Erikss orTet al.lll9 83) in temperature 
until the emission parts are removed and the line cores are fitted. 
This leads to a plateau temperature of T v = 4, 500 K, instead 
the value of T v = 5,500 K assumed in the original model. To 
maintain the physical conditions of the chromosphere, we also 
scaled the electronic number density in order to keep the elec- 
tronic pressure constant. Mor eover, we also i ncreas ed the mi- 
croturbulent velocity profile of E riksson et al.l dl983)'s chromo- 
sphere by 20 % in order to enlarge the gaussian core of the lines. 
The hybrid model with T p = 4, 500 K a nd larger microturbul ent 
velocity is compared with the model of Eriksson et al.l d 1983b in 
Fig.0 As shown in Fig. [12] we satisfactorily reproduce the line 
cores, but only when applying a macroturbulent parameter of 4 
and 1 1 for the Ca n/849 nm and Ca n/854 nm lines, respectively 
to the three models. However, we note that our model does not 
perfectly reproduce the line core for the Ca 11/866 nm line. For 
strong lines such as the Can/ triplet, the weak projected rota- 
tional velocity of /3 Cet has no effect on the line cores. We note 
the difference in the line core between an interpolated MARCS 
model without a chromosphere and the hybrid model with a 
chromosphere. This difference can be used as an indicator of a 



available at: http://masics.basn.obs-mip.fr 



chromosphere an d stellar activity (e.g. lAndretta et al. 200l and 
Busa etalJl2007h . 



6. Conclusion 

By analyzing interferometric measurements derived from 
CHARA/VEGA data, we have determined the physical extents 
of the chromosphere of late-type giant stars. For the first time, we 
have measured the position of the area in the atmosphere where 
Can triplet and H a lines formed and found that these chromo- 
spheres are relatively extended, between 16% and 47% of the 
stellar radius of the corresponding photospheres. Surprisingly, 
we have found that four of the program stars exhibit anoma- 
lous phases of the interferometric fringes that can be interpreted 
as an inhomogeneity in the chromospheric part of the atmo- 
sphere. More observations with a more complete (m, v) coverage 
are needed to explain the detected asymmetry of these chromo- 
spheres. 

In addition to this chromospheric study, we have derived 
the fundamental parameters of seven cool giant stars. From the 
stellar radius measured in the continuum and the corresponding 
bolometric flux, we have deduced the effective temperatures of 
the program stars, which are in excellent agreement with T e g 
scale laws found in the literature for such K giant stars. 

We have also developed a semi-empirical model atmosphere 
of (3 Cet. This model has helped us to define the limb-darkening 
law to apply to the core of the Ca 11 triplet lines. We have found 
that it is similar to the one found for the Sun. We have used 
this model to define where the cores of the Ca 11 triplet lines are 
formed by analyzing the contribution functions. We have been 
able to confirm that these lines are formed at the mean chro- 
mosphere temperature and that the Ca n/849 nm line is formed 
significantly deeper within the atmosphere than the two other 
Can triplet lines. This result is in very good agreement with 
the interferometric measurements for /3 and 77 Cet. Finally, we 
have computed synthetic Can line profiles, which have been 
compared to an observed high resolution spectrum of /3 Cet. 
We have found that the temperature plateau (as proposed by 
lEriksson et al.|[l983l) of the chromosphere model is too high to 
reproduce the line cores. 

On the basis of these first interferometric results on the spa- 
tial structure of the chromosphere of K giant stars, it now seems 
possible to construct reliable chromospheric models by combin- 
ing spectroscopic and interferometric observations. 
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Fig. 9: Electronic temperature, electronic number density and 
microturbulent velocity profiles as a function of the integrated 
mass density log m of f3 Cet used to reproduce Ca n triplet line. 
Dashed line is for model of Erik sson et al.l d!983l) . Full line is 
for a hybrid model betwe en an interpolated MA RCS model and 
a scaled chromosphere of Erik sson et alJ dl983l) 's model. 



Acknowledgements. VEGA is a collaboration between CHARA and OCA / 
LAOG / CRAL / LESIA that has been supported by the French programs PNPS 
and ASHRA, by INSU and by the Region PACA. The CHARA Array is operated 
with support from the National Science Foundation through grant AST-0908253, 
the W. M. Keck Foundation, the NASA Exoplanet Science Institute, and from 
Georgia State University. This research has made use of the Jean- Marie Mariotti 
Center LITpro and SearchCal services co-developed by CRAL, LAOG and 
FIZEAU, and of CDS Astronomical Databases SIMBAD and VIZIER. Finally, 
we thank the anonymous referee for providing pertinent advices that helped im- 
proving this paper. 



References 

Andretta, V, Busa, I., Gomez, M. T, & Terranegra, L. 2005, A&A, 430, 669 

Ayres, T. 2010, Mem. Soc. Astron. Italiana, 81, 553 

Baines, E. K., Dollinger, M. P., Cusano, F, et al. 2010, ApJ, 710, 1365 

Berio, P., Stee, P., Vakili, F, et al. 1999, A&A, 345, 203 

Bonneau, D., Clausse, J., Delfosse, X., et al. 2006, A&A, 456, 789 

Brown, A. & Carpenter, K. G. 1984, ApJ, 287, L43 



1 .0 
0.8 
0.6 
0.4 
0.2 



0.0 



Call 849.8 nm 
Call 854.2 nm 
Call 866.2 nm 
Call continuum 



0.0 



0.2 



0.4 0.6 

sin 9 



0.8 



Fig. 10: NLTE computed limb-darkening laws for Can triplet 
line cores (A/1 = 0.25 A) and continuum with the hybrid model 
for j3 Cet as a function of sin 9 = - ji 2 . 



Busa, I., Aznar Cuadrado, R., Terranegra, L., Andretta, V, & Gomez, M. T. 2007, 
A&A, 466, 1089 

Carlsson, M. 1986, Uppsala Astronomical Observatory Reports, 33 
Carpenter, K. G., Brown, A., & Stencel, R. E. 1985, ApJ, 289, 676 
Chen, B., Vergely, J. L., Valette, B., & Carrara, G. 1998, A&A, 336, 137 
Claret, A., Diaz-Cordoves, J., & Gimenez, A. 1995, A&AS, 114, 247 
Cuntz, M. 1990a, ApJ, 349, 141 
Cuntz, M. 1990b, ApJ, 353, 255 

di Benedetto, G. P. & Rabbia, Y. 1987, A&A, 188, 114 
Diaz-Cordoves, J., Claret, A., & Gimenez, A. 1995, A&AS, 110, 329 
Dyck, H. M., Benson, J. A., van Belle, G. T., & Ridgway, S. T. 1996, AJ, 111, 
1705 

Eaton, J. A. 1993, ApJ, 404, 305 

Eriksson, K., Linsky, J. L., & Simon, T. 1983, ApJ, 272, 665 
Gustafsson, B., Edvardsson, B., Eriksson, K„ et al. 2008, A&A, 486, 951 
Flaisch, B. M. 1987, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 

291, Cool Stars, Stellar Systems and the Sun, ed. J. L. Linsky & R. E. Stencel, 

269 

Hall, J. C. 2008, Living Reviews in Solar Physics, 5, 2 
Houdashelt, M. L., Bell, R. A., & Sweigart, A. V. 2000, AJ, 1 19, 1448 
Judge, P. G. 1987, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 291, 
Cool Stars, Stellar Systems and the Sun, ed. J. L. Linsky & R. E. Stencel, 294 
Judge, P. G, Jordan, C, & Rowan-Robinson, M. 1987, MNRAS, 224, 93 
Kovacs, N. 1983, A&A, 120, 21 
Linsky, J. L. & Haisch, B. M. 1979, ApJ, 229, L27 
Luck, R. E. & Challener, S. L. 1995, AJ, 1 10, 2968 

Massarotti, A., Latham, D. W., Stefanik, R. P., & Fogel, J. 2008, AJ, 135, 209 
McWilliam, A. 1990, ApJS, 74, 1075 

Meilland, A., Stee, P., Vannier, M., et al. 2007, A&A, 464, 59 
Merle, T., Thevenin, F, Pichon, B., & Bigot, L. 201 1, ArXiv e-prints 
Montes, D., Fernandez-Figueroa, M. J., De Castro, E., et al. 2000, A&AS, 146, 
103 

Morel, P. 1997, A&AS, 124, 597 

Mourard, D., Clausse, J. M., Marcotto, A., et al. 2009, A&A, 508, 1073 
Mozurkewich, D., Armstrong, J. T., Hindsley, R. B., et al. 2003, AJ, 126, 2502 
Perryman, M. A. C, Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49 
Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168 
Reimers, D. 1982, A&A, 107, 292 

Richichi, A., Percheron, I., & Davis, J. 2009, MNRAS, 399, 399 
Simon, T., Linsky, J. L., & Stencel, R. E. 1982, ApJ, 257, 225 
Skrutskie, M. F, Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 
ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2005, ApJ, 628, 
453 

van Belle, G. T, Lane, B. F, Thompson, R. R., et al. 1999, AJ, 1 17, 521 
van Leeuwen, F. 2007, A&A, 474, 653 
White, O. R. 1962, ApJS, 7, 333 



10 P. Berio et al.: Chromosphere of K giant stars 




■4 -3-2-10 1 2 -5 -4 -3-2-10 1 2 -5 -4 -3-2-10 1 2 

log rn log rn log m 



Fig. 1 1 : Contribution function s to both the co r e (thin lines) and the continuum (thick lines) of the Ca n triplet line profiles in (i Cet. 
The dashed lines represent the Eri ksson et al.l (119831) model atmosphere. The full lines are for the hybrid model. 




Fig. 12: Fits of the NLTE Can triplet line pr ofiles in /3 Cet . The d otted line is the observations from NARVAL at TBL with a 
resolution of ~ 65, 000. The dashed line is the lEriksson et al.l {1983) model atmosphere with a chromosphere (T p = 5, 500 K). The 
thin full line is the interpolated MARCS model. The thick full line is our hybrid model with a chromosphere (T p = 4, 500). The grey 
lines below the spectral line profiles show the rotation profiles with a projected rotational velocity of 5 km s _1 . 



